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convincingly  demonstrates  the  potential  of  the 
consistent  reference  state. 


moments  of  a molecular  system  may  be  evaluated  directly  by 
the  Polarization  Propagator  II].  Spectroscopic  properties 

guadrupole  moments,  radiative  lifetimes,  Verdet  constants, 
and  frequency  dependent  polarizability  are  easily  understood 
and  calculated  within  the  propagator  formalism.  Propagators, 
or  the  Green's  functions  as  response  function  to  an  external 
field  probe,  have  originally  found  successful  applications 
in  field  theory  (2],  and  statistical  mechanics  [3].  The 

function  referring  to  directly  measurable  quantities  is 
particularly  appealing  for  all  areas  of  spectroscopic 
studies  and  has  been  implemented  into  statistical  physics  by 
Zubarev  [4]  and  into  quantum  chemistry  by  Linderberg  and 
Ohrn  (1).  Calculations  of  excitation  energy  by  this  approach 

states.  Numerical  inaccuracy  due  to  taking  the  difference  of 
two  very  large  numbers  (state  energies)  to  obtain  excitation 

optimized  reference  state  to  a polarization  propagator  must 
form  a basis  for  a balanced  description  of  both  ground  and 
excited  states  in  order  to  produce  accurate  excitation 
energies,  and  their  associated  properties.  Early  quantum 
chemical  calculations  of  the  Polarization  Propagator  have 


applied  the  linearized  tine-dependent  Hartree  Pock  (TDHP) 
[S]  or  the  Random  Phase  Approximation  { RPA ) , [6],  and  the 
Tamm  Dancoff  Approximation  ( TDA ) , [7],  using  an  optimized 

single  determinantal  (Hartree  Pock)  reference  state  (8).  Por 
a complete  set  of  particle-hole  (a*aw)  and  the  hole-particle 
(ajla  1 operators,  a unique  set  of  excitation  operators  for 
the  RPA  has  been  given  by  Dalgaard  (91, 


(1) 


qJJ  (Reference  state?  - |Excited  state,  k>  (2) 

where  q^*  is  a linear  combination  of  particle  hole 
operators,  and  qk  its  adjoint,  and  are  the 
complex  expansion  coefficients.  For  a consistent 
approximation  for  RPA  one  must  then  also  fulfill  the 
condition 


| Reference  state?  - 0 (3) 

The  failure  of  the  Hartree  Pock  state  to  satisfy  the 
above  condition,  has  led  to  several  approximations  aimed  at 
remedying  this  inconsistency.  A higher  RPA  which  uses 
multideterminantal  TDHP  by  constructing  a linear  combination 
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Antisymmetrized  Geminal 


(GAGP)  wavef unction . 


is  of  same  form  as  the  N particle  projected  Bardeen  Cooper 
Schrieffer  (BCS)  function  used  in  the  theory  of 
superconductivity  [16|  and  nuclear  physics  [17J.  While  the 
GAGP  in  practice  has  not  proven  to  be  an  exact  consistent 
reference  to  the  Polarization  Propagator  in  the  RPA  the  GAGP 
is  a better  consistent  reference  state  to  the  Generalized 
Tamm  Dancoff  Approximation  (GTDA),  than  is  the  single 
determinantal  SCP.  The  GAGP  based  GTDA  permits  us  to 
identify  the  eigenstates  obtained,  as  spectroscopic  excited 
states.  Such  an  approximation  amounts  to  setting  y - 0 in 
equation  (1),  and  it  has  been  conjectured  that  the 

smaller  than  in  an  uncorrelated  reference  state  as  HP.  To 
test  the  validity  of  such  a viewpoint,  comparative  studies 
were  carried  out  for  Be,  LiH,  and  Lij.  These  results  give  a 
strong  indication  that  indeed,  the  RPA  results  are  much 
closer  to  the  TDA  when  the  reference  state  is  an  AGP  rather 
than  a HP  state.  These  results  are  tabulated  in  Chapter 
Pour.  These  results  further  indicate  that  the  triplet 
instability  problem  is  greatly  improved  by  using  AGP  over 
HP.  The  Generalized  form  of  the  AGP  (GAGP)  contains  a single 
determinantal  independent  particle  function  that  could 
accommodate  core  electrons  which  do  not  require  detailed 
correlation  or  accommodate  any  unpaired  electron  (for 
doublet  spin  systems)  which  cannot  be  included  in  the 
geminal  of  the  AGP.  The  use  of  the  generalized  form  often 


proves  to  be  more  efficient  in  comparison  with  the  fully 
correlated  AGP.  There  is  usually  no  loss  in  quality  of 
results  (by  freezing  the  core  orbitals  to  configurational 
variation)  (18).  The  case  of  the  open  shell  systems  is 
complicated  by  the  fact  that  even  though  the  adjoints  of  the 
excitation  operators  satisfy  the  vacuum  condition,  the  GTDA 
fails  to  be  a consistent  approximation  for  the  open  shell 
(spin  doublet  ground  state)  GAGP  reference  state.  This  has 
been  elucidated  in  Chapter  Three  and  the  problem  arises 
primarily  due  to  the  spin  degeneracy  and  the  need  to  span 
the  full  spin  space.  In  this  case  one  can  only  expect  to 
obtain  numerical  results  using  the  full  RPA  expression.  In  a 
paper  by  SwanstriSm  and  Jorgensen  (191,  the  authors  have 
evaluated  excitation  energies  of  open  shell  systems  using 
the  perturbative  propagator  formalism  where  they  have 
treated  the  excited  state  with  an  ct  unpaired  electron  as 
their  zeroeth  order  wavefunction  and  the  configuration  with 
an  unpaired  electron  with  a 6 spin  as  a second  order 
correction  to  the  wavefunction.  with  the  same  strategy,  one 
may,  in  the  first  approximation,  simply  consider  the 
configuration  of  excited  states  to  be  consisting  of  an 


unpaired  electron  with  a spin  only  (which  is  then  correct  up 
to  first  order  in  the  perturbation  expansion).  It  is  an 

this  approximation.  We  have  programmed  the  GTDA  for  this 
approximation  and  used  it  to  evaluate  excited  states  of  the 
Li  atom.  The  first  few  excitation  energies  are  surprisingly 


close  to  experimental  values,  but  one  needs  to  check  such  a 
procedure  tor  larger  systems  before  any  conclusions  can  be 
drawn.  It  nonetheless  provides  an  inspiring  first  step 
towards  developing  a generally  valid  theory,  although  a more 
involved  numerical  procedure  is  outlined  in  Chapter  Three 
together  with  results  of  these  preliminary  calculations. 


CHAPTER  ONE 

POLARIZATION  PROPAGATOR 


Background 

We  nay  describe  by  a double-time  Green's  function  [1) 
the  causal  probability  amplitude  of  the  density  fluctuation 
function  p(t),  of  a pure  state,  as  it  responds  to  a weak 
external  field.  The  double-time  Green's  function  is 
expressed  as, 

G(t,t')  • <<p(t),p(t')» 

- i8(t-t')<0|p+(t)p(t')|0> 

- i6!t'-t)<0|p+(t’)p(t)|0>  (1-1) 

where  9 is  the  Heaviside  step  function,  such  that 


and  |0>  is  the  ground  state. 

Invoking  the  resolution  of  the  identity  in  equation  (1-1), 
with  a complete  set  of  discrete  eigenstates  |m>,  of  the 
hamiltonian  H (equation  (1-4)),  permits  us  to  write,  (1), 

«p(t)jp(t')»  - l [ie(t-t')<0|p(t)|mxm|p(t')|0> 


(1-2) 
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spin  orbital  i,  and  Cj,  Cj  are  coefficients  for  the 
excitation  operators.  The  upper  signs  corresponds  to  an 
excitation  operator  for  a spin  singlet  state  and  the  lower 
signs  to  an  operator  for  a triplet  state.  The  most  general 
one  particle  excitation  operator  may  then  be  constructed  by 
a linear  superposition  of  the  q+  ' s , and  the  g 's  [9].  The 
new  set  of  excitation  operators  (Qv+,  Qv)  may  be  written  as, 


[ »*  « ] ■ [ «*  « ] | ■ > 

(1-211 

with  (Z,Y)  as  expansion  coefficients  for  general  excitation 
operators. 

For  { Q ,Q  } to  form  a set  of  excitation  operators, 
Dalgaard  (91  proved  the  following  theorem.  Let  |0>  and  |X> 
be  arbitrary  N-electron  states,  which  satisfy  <0|+n>  * 0 and 
<0|X>  - 0.  If  the  matrix  <(g,g'l'l>  is  non  singular,  there  is 


5 ( <3 


(1-22) 


|X>  - q£|o> 

Qx|0>  - 0 


(1-23) 


In  the  above  theorem  l+g)  is  an  N-electron  state  in  the  Fock 
space  and  can  be  written  as 


IV 


(1-24) 


The  polarisation  propagator  F{ E)  is  formally  expressed  as 


Incorporation  of  the  propagator  expression  of  equation  (1- 
18)  into  equation  (1-25)  yields  the  following  matrix 
equation, 


with  the  notation, 


A - <|(q.!H,q+]]|> 

B - <|lq,[H,q]]|> 

and  the  set  |q+,q)  satisfy  the  relations. 


(1-27) 


<|l9  . q+)|>  - X - - <|(q+,  qj | > 

<|lq+.  q+l|>  - o - < I Iq  , qj|> 

Diagonalization  of  propagator  matrix  is  accomplished  by  a 
series  of  transformations  [14,28],  so  one  may  express  the 
propagator  as. 


The  eigenvalues  a>,  of  the  propagator  are  solutions  of  the 
eigenvalue  equation, 


Even  though  the  propagator  matrix  on  the  left  hand  side  of 
equation  (1-30)  is  itself  hermitian  and  has  real  eigenvalues 

[29] ,  this  generalized  eigenvalue  equation  may  have  complex 

eigenvalues.  This  is  related  to  the  stability  of  the 
optimized  reference  state  with  respect  to  symmetry  breaking, 
which  will  be  discussed  in  Chapter  Two.  If  the  reference 
state  is  stationary  with  respect  to  the  variational 

parameters  then  the  propagator  matrix  is  in  fact  the  Hessian 

[30]  of  the  variational  problem  (see  Chapter  Two). 
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Ideally  we  require  the  operators  [Q+J  to  be  such  that  their 
adjoints  {0}  annihilate  the  reference  state. 


Qm|0> 


(1-35) 


/er  well  known  that  the  Hartree  Pock  reference 
not  satisfy  this  vacuum  condition  (1-35)  and  is 
consistent  reference  state  to  the  RFA  (10,14,15). 


only  if  Y - 0 in  equation  (1-34).  One  may  set  Y-0  as  an 
approximation  to  the  RPA,  known  in  the  literature  as  the 
Tamm  Dancoff  approximation  ( TDA ) . 


The  inconsistency  of  the  RPA  is  rooted  in  the 
interaction  of  the  reference  state  with  excited  states 
through  the  hamiltonian  (which  is  symbolized  by  the  B matrix 
in  the  propagator).  Consistency  requires  that  the  B matrix 
be  zero  identically.  In  the  absence  of  any  reference 
function  which  fulfils  this  requirement  exactly,  one  may  as 
an  approximation  set  the  B blocks  to  zero  and  obtain  the 
Tamm  Dancoff  Approximation  as  a 
the  RPA.  Setting  the  B block 


consistent  approximation 


eliminated  and  the  excitation  operators  for 
then  simple  particle-hole  operators. 


Where  N is  the  total  number  of  doubly  occupied  orbitals  and 
2S  is  the  total  number  of  basis  functions.  It  is  obvious 
that  the  Q^(TDA)  satisfy  the  vacuum  condition  since  the  i 

annihilate  the  reference  state. 

With  this  background,  we  will  now  look  into  a more  general 
and  correlated  wavefunction  which  may  be  a suitable 
reference  function  for  the  polarization  propagator.  As 
mentioned  in  the  introduction,  the  Generalized 
Antisymmetrized  Geminal  Power  wavefunction(GAGP)  is  such  a 
reference  state  [14,15].  While  in  the  past  few  years  the 
AGP-TDA  ansatz  has  produced  some  interesting  results  [32], 
in  this  project  we  have  developed  and  investigated  several 
aspects  of  the  AGP.  We  have  developed  the  generalized 


also  made 


comparative  studies  of  the  GTDA  and 
measure  of  consistency  of  this  ansatz. 
the  AGP  ansatz  are  compared  with  t[ 
reference  state.  The  next  chapter  coni 


PA  to  estimate  the 
calculations  with 


a discussion  of 


Generalized  Antisvmnetrized  Geminal  Power  Wavefunction 


The  AGP  function,  formulated  as  an  N-particle  projected 
BCS  (Bardeen,  Cooper,  Schrleffer)  function  and  alternatively 
known  as  the  Scaff roth-condensed  pair  wavefunction  [17],  has 
been  applied  successfully  in  the  theory  of  superconductivity. 
Efforts  directed  towards  obtaining  a proper  reference  state 
for  the  polarization  propagator  have  led  to  suggestions  by 
tinderberg  and  Ohrn  [14],  and  Ohrn  and  Linderberg  [15]  for 
employing  AGP  as  the  reference  state  for  polarization 
propagator.  A generalized  AGP  function  was  also  established 
by  Goscinski  and  Weiner  [33]  to  be  the  correct  reference 
state  for  the  self  consistent  one  particle-hole  propagator 
( SC1PHP ] . Furthermore,  because  of  its  sparse  and  well 
structured  first  and  second  order  density  matrices,  the  AGP 
function  has  been  explored  by  Coleman  [34],  Nakamura  [35], 
and  others  in  connection  with  the  N-representability 
problem.  Certain  early  attempts  with  the  AGP  reference  state 
had  been  rather  discouraging  because  of  inefficient 
computational  schemes.  However,  in  recent  times 
sophisticated  optimizing  strategy  built  on  exploiting  the 
particularly  well-structured  characteristics  of  first  and 
second  order  density  matrices  of  the  AGP  has  proved  the 
wavefunction  gives  an  efficient  description  of  electron 


correlation. 


Functional 


The  basic  information  in  an  AGP  wavefunction  is 
contained  in  a single  geminal  G(ij),  which  represents  a pair 
of  electrons  i,j.  The  geminal  may  be  expressed  as  a 
superposition  of  all  two  electron  determinants 
weighted  by  factors  g^.  He  may  express  G+,  creating  a 
geminal  on  the  vacuum  as, 

0*11)1  - ! «„{  - •Jli)-Jl)l  } U-n 


where  2s  is  the  total  number  of  functions  in  the  basis  set. 

and  are  natural  spinorbitals  (nso's)  that  contain  the 
same  spatial  function  but  with  a and  0 spins  respectively. 
Thus  for  a system  with  N(even)  electrons,  an  AGP  function 
is  an  antisymmetrized  N/2  power  of  the  geminals  (14,15), 

AGP  - (Cn/2I_1(G+)N/2  lV»C>  <2"2> 

where  G(i,j)  may  also  be  written  in  the  natural  spinorbital 
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[36], 


also  considered  two  different  geminals  in  the  AGP  with 
interesting  results  for  the  ground  state  of  Be.  In  this  case 
the  total  energy  of  the  the  system  is  the  same  as  the  MCSCF 


total  energy  with  the  same  basis.  But  it  is  important  to 
note  that  a true  AGP  function  is  a power  function  of  a 
single  geminal.  It  is  often  convenient  to  attach  to  the  AGP 
function,  an  independent  particle  factor  (IPF)  ♦ , that  can 
account  for  molecular  orbitals  which  are  essentially 
uncorrelated,  such  as  the  core  and  some  inner  shell 
electrons  in  a reasonably  large  system.  The  IPF  component 
also  becomes  essential  for  calculations  involving  open  shell 
systems  (with  odd  number  of  electrons)  to  accomodate  the 
unpaired  electrons.  The  Generalized  Antisymmetrized  Geminal 
Power  ( GAGP ) wavefunction  for  N electrons  is  formulated  as 
an  antisymmetrized  product  of  AGP  and  IPF  components, 
containing  k geminals  and  N-2k  IPF  functions. 


GAGP  - (c  -UN  - 2k)  l)_1  0AS  { G*  }k  ajk+1 . . . . aj | Vac> 

(2-7) 

The  AGP  wavefunction  contains  a very  large  number  of  config- 
urations, (jj)  with  much  smaller  number  ( s ) of  the  variatio- 
nal parameters  g.  The  energy  fuctional  E(C,g),  which  is  a 
function  of  a set  of  orbital  coefficients  C,  and  a set  of 
geminal  expansion  coefficients  g,  is  represented  by 


MSP  I H | GAGI 

CGAGP  [GAGP> 


(2-8) 


Since  the  hamiltonian  consists  of  at  most  two-particle 
interaction  terms,  one  can  express  the  above  energy 
expression,  by  means  of  a reduced  hamiltonian  K„,  and  second 
order  density  matrix  ^ (381,  in  the  following  form. 


Equation  (2-9)  is  much  more  attractive  than  equation  (2-8) 
computationally  because  of  reduced  dimension,  however  (2-9) 
cannot  be  used  unless  the  second  order  density  matrix 
corresponds  to  states  in  N-dimensional  Hilbert  space.  This 
is  known  in  the  literature  as  the  N-representability  problem 

and  we  can  use  equation  (2-9)  to  compute  the  energy  quite 
simply.  We  use  second  quantized  operators  to  describe  the 
reduced  hamiltonian  and  2-density  as  follows. 

The  reduced  hamiltonian  K„  has  the  form. 


,> ) 


(2-9) 


is  the  second  order  density  matrix  of  the  GAGP. 


<ij||kl>)  ajaja^ 


(2-1 


The  only  possible  types  o£  non-zero  second  order  density 
matrix  elements  [34]  are  summarized  below  : 


( 2 ) + , , SU.J> 

Da  • <GAGP | a ^a-a-a j |GAGP>  - g^  aN-l 

■ / S|J 

1 J 

is  a square  matrix  of  order  ( 2 )+  ^ 

d'2>  . <GAGP | a jB^a^a j |gagp> 

<GAGP I a ja Ja ja J | GAGP> 

<GAGP  |a taJaj-a^  |GAGP> 

<GAGP  | aiataja^ |GAGP> 


is  a diagonal  matrix  of  order  (jj) 

D^2’  ■ rGAGPlataJa^jlGAGPi 
<GAGP|a^a£aBa5|GAGP> 

<GAGP|a*a|aBaj|GAGP> 

<GAGP|a^a*abaj|GAGP> 


(2-13) 


“(N-2k)  order  diagonal 


■Ja‘ 


■ <GAGP|»*a;»bllc|GAGP> 
<GAGP | a J«5«Bac I SAGP> 
<GAGP|,*.gaB,c|GAGP> 
<GAGP|a|aJ,b.j|GAGP> 


are  ordered  such  that  l<i,j<k,  k+l<b,c<(N-2k) . 


E(C,  9)  - ««,  * *0>  D<2>+  D<2>+  Kd  . 1 

■si  [('■.'«  * k2\, 

* i <V,,  *j'V„] 


E(X,g)  - 

"M'” • *:■,>)• 


(2-16) 


,Xcs>  is 


E(X,g)  - E(0,g)  < 


V.7X  E(0,g)  + i Xt.V^E(0,g)  X (2-17) 


; V2E(X,g)  > 0.  we  thus  obtain  the  r 


*n<w«  - “Ki- 

When  the  commutator  in  (2-19)  is  expand 


5kl(\»in'9>  " Re  Tr|  laIeiHKH.D<2>]|  (2-20) 

vkl(Xmin'9»  - Re  Trf  I l^k^r><^r^1l^RH'^^>^  (2_21> 

”ki'W®»  - Re  Tt<  :<Vti>KH'Dl2,'iw  <2-22> 

Tkl‘x»in'9>  ■ Rb  ( 5IRH-D(2,'krlr) 

- Re  r5,s,t["CH)K«t  “stir  - “krst'Vstlr] 


(2-23) 


”K»  - *w-x  - •>,'  * wi\  * - 

■j  [ v*;*q  - •;*»'  * “■  v'X  * ■;v)]vo 

■»  i * »!,(«*  w>.  - *:v  * **■  - 

- •>,>  * * -xor  •••>”’) 


*«..  ■ *X  ' ‘5*5 

- ■;*»  * "i*s  (2-2H 


*■(?’-  *£,(“  V-S,  - w * “■  W*U  * v)  - 
-i  [ ,!„[«■  v*;,  - v * *■■  v-;,  * 

>«  (■  * *!„(••  * >“  >■<,>*;.  * *,.0  - 

(2-26) 


V2E<X,g)ij|(1|X-*0  " 

rr[{(  X*,  - x^C  - xkl)  + ( xjj  - xkl)xB(  x^  - xkj) 

"'I  ((  Xtj  * *Il  * XklK  + < Xlj  - Xij^  Xkl  - Xkl)KH 
+ *H(  xtj  - xji)(  xIi  - hJ  * kh(  xIi  - xici)(  x+a  - xij)))D<2>] 


Introducing  the  following  notation. 


Aijkl  - Tr[([xJj  .[kh  - Jtjj]]  + [x+j  ,[kh  , x+3]])d(2>] 

Bi  jkl  - «t([*I3  ,[kh  , Xkl]]  ♦ [xj,  ,[kh  , Xtj]])D<2>] 

(2-28) 


Equation  (2-27)  may  be  written  in  symbolic  form  as 


} ijkl  | X-»0  ' Z ( Aijkl  Bi  jkl  " 

- Re(A  - B) j jkl 


Ai jkl) 


in  equation  (2-28)  we  have  used  the  fact  that  the 
expectation  value  of  the  adjoint  of  an  operator  is  the 
complex  conjugate  of  the  expectation  value  of  the  operator. 
2)Variation  with  respect  to  Im  X. j and  Im  Xkl  gives, 

72E(X,g)ijki|x<0  - 

Tt  [”  Xij  + Xi  j^XH  ^ Xkl  + Xkl^  + ^ Xkl  + Xkl^XH^  Xij  + Xij^ 
"-J  ((  Xij  + Xij)(  Xkl  + XklK  + ( Xij  + Xij^  <1  * Xkl)KH 
+ kb(  Xtj  ♦ x13)(  *+,  ♦ xlk)  ♦ KH(  X£l  ♦ xkl)(x*j  * x.3))}d<2»] 

(2-30) 


Equation  (2-2 


written  symbolically 


Aijkl) 


3)  Variations  with  respect  to  Re  Xjj  and  Im  Xkj  yields  the 
tollwing  expression 


v2E(X.g)lj|(l|x^o  - 

»[*  <<  <t,  - *.)>*„(  <1  * »„)  *<*;>*  >«.)■,(  - <«) 

'2  Xij  * Xij^  Xkl  + XkPXH  + ^ Xij  ■ XiJ^  Xkl  + Xkl^KH 
♦ kh(  *I3  - x4j)  ( xj,  + Xlk)  ♦ XH(  ^ ♦ xkl)  ( X*,  - Xij))}p<2>] 
(2-32) 

Which  has  the  following  symbolic  represention 
V E( X, g) j jkl  |k_,0  - j ( Aljkl  - B1Jkl  + Bjjkl  - AjJkl) 


4)  We  obtain  the  hessian  with  respect  to  1m  X_  and  Re  Xkk 


”[*{< 
"-!(( 
* >,<  ■ 


!,ijkl|x-»0  " 
Xt,  * ’‘ijM 

■»**»>< -I 


'■■I1 ' ■ 'ki  -ti)",'  xij  - *ij) 

l)«.*  C*I3 

| (4j  ♦"!)%'=’] 


which  may  rewritten  in  symbolic  form  as 


2 ( Ai jkl  " Bi jkl 

lm(A  - B)i:jkl 


(2-35) 


Thus  we  find  that  a complete  singlet  hessian  consists  of 
matrix  elements  of  second  order  variations  with  respect  to 
real  components  of  X,  and  imaginary  components  of  X, 
including  complex  variations  which  contains  mixed 
components.  Our  optimization  procedure  at  the  moment  makes 
use  of  only  real  variations  of  the  singlet  block  and  utilize 
only  equation  (2-26)  to  obtain  the  hessian.  Explicit 
expression  for  the  hessian  elements  are  included  in  appendix 

The  complete  hessian  is  blocked  according  to  each  type  of 
spin  symmetry  variations  such  as  singlet,  triplet,  etc. 
However,  since  we  want  to  retain  the  singlet  of  the  ground 
state,  we  do  not  use  other  blocks  of  the  hessian.  This 
restriction  sometimes  leads  to  problems  regarding  RPA 
calculations  where  instabilities  occur  in  generating  triplet 
excitation  energies.  Zn  the  same  way  one  can  block  the 
hessian  according  to  spatial  symmetries  such  as  , say,  E,  0 
etc.  for  linear  molecules,  and  S,  P etc.  for  atoms.  A 
variation  of  the  sigma  type  retains  the  original  symmetry  of 
the  molecule  in  its  reference  state,  which  is  usually  a 
sigma  state.  Other  types  of  variations  will  produce  an 
optimized  state  with  its  symmetry  corresponding  to  the  cross 


product  of  the  symmetry  of  the  infinitesimal  generators  and 
that  of  the  reference  state.  For  example,  generators  of  n 
symmetry  acting  on  a sigma  reference  state  would  produce  a 
n state.  Since  we  only  optimize  our  reference  with  respect 
to  the  sigma  generators,  'instabilities'  in  non-sigma  blocks 
of  the  propagator  matrix  may  occur  in  calculations 
indicating  that  symmetry  broken  solutions  are  possible  with 
lowering  of  energy. 

The  Excited  States  Of  The  GAGP 

The  excitation  energies  and  the  excited  states  of  the 
GAGP  are  obtained  by  the  excitation  operators  as  discussed 
before,  and  a particular  set  of  N-particle  operators  known 
as  the  diagonal  operators.  The  basis  set  (q+,q)  of  one 
particle  excitation  operators  may  be  defined  to  be  the 
following, 


(2-36) 


nomalization  constant 


(2-37) 


9i  > 


(2-38) 


36 


Nj  - <0|atai  + afaj|0> 

The  N-particle  diagonal  operators  are  essential  lor  totally 
symmetric  singlet  excitations.  These  operators  are  of  the 


(2-39) 


The  superscripts  s,t  in  (2-35)  refer  to  singlet  and  triplet 
symmetries.  Since  one  may  not  excite  into  an  IPF(occupied) 
spinorbital  the  indices  i and  j are  such  that  1 $ l £ k,  and 

5 i S N*  The  adjoints  of  operators  annihilate  the  reference 


qj^’lGMSP)  . o 


r all  possible  i,j 


While  in  practice,  an  exact  propagator  calculation  using  a 
complete  set  of  operators  is  unfeasible,  certain  significant 
properties  inherent  in  the  optimized  GAGP  function  allow 
meaningful  interpretation  of  the  propagator  with  the  above 
truncated  set  of  operators.  Weiner  [431  has  proved  that, 
when  a reference  state  satisfies  the  following  relations. 


:|[H,  X+Y) |Re£> 


(2-41) 


<Ref|(H,  X] |Ref>  - 0 (2-42) 

Por  all  operators  X,'  Y belonging  to  an  operator  subset  (V) , 
it  is  possible  to  construct  a projected! model)  hamiltonian 
superoperator,  PHP,  which  acting  on  the  entire  hilbert  space 
of  operators,  leaves  the  subspace  spanned  by  the  basis  set 
of  V invariant.  In  such  a case  the  eigenvectors  and 
eigenvalues  obtained  by  the  action  of  PHP  and  the  full 
hamiltonian  on  all  such  operators  X,  Y are  identical.  If  the 
reference  state  is  a GAGP,  the  projection  operator  is 


P - ^ j k x 1ij|GAGP><GAGP|qq+|GAGP>'1<GAGP|gkl 

+ |GAGPXGAGP  | (2-43) 

The  GAGP  function  does  not  satisfy  the  relations, 

< GAGP | (H,  q^qj^ilGAGP)  - 0 

<GAGP | [ H,  qjjqJj! |GAGP>  - 0 

<GAGP  1 1 H,  qi;jqkl]|GAGP>  - 0 (2-44) 

However,  the  optimized  GAGP  function  does  satisfy  the 
generalized  Brillouin  condition, 


<GAGP| (H,q+) |GAGP> 


'|[H,  q] |GAGP> 


(2-45) 


and  a very  significant  relation  that  was  proved  by  Jensen, 
Weiner  and  Ohrn  127],  and  Sangfelt,  Goscinski,  Kurtz  and 
Blander  (44],  is, 

<GAGP  Iq^  jqj^H  |GAGP>  - *i|c*j1<GAGP  |H  |GAGP>  (2-46) 

A consequence  of  the  GAGP  fulfilling  relations  (2-45),  (2- 
46),  is  that  if  the  following  propagators  are  constructed, 
based  on  the  full  hamiltonian  (2-47),  and  the  model 
hamiltonian  (2-48), 


I<GAGP | [q, [H,q+] ] |GAGP> 
<GAGP|(qt(H,q+]] | GAGP > 

r<GAGP| (q, [ PHP,q+] ] |GAGP> 
I < GAGP | Iqt( PHP , q+ ] ] |GAGP> 


<GAGP | tq, (H,q) ] | GAGP >1 
<GAGP|(qt[H,q]] |GAGP>J 

<GAGP | [q, [ PHP,q] ] |GAGP>1 

<GAGP I Iqil PBP,q) ] | GAGP> I 
(2-48) 


the  diagonal  blocks  of  G and  Gr  are  identical,  so  that  the 
generalized  TDA  based  on  a model  hamiltonian  is  the  same  as 
that  based  on  the  full  hamiltonian.  The  subset  of  operators 
used  to  construct  the  GTDA  fulfills  conditions  expressed  by 
(2-41,2-42).  Consequently,  it  follows  that  if  we  consider 
the  excitation  operators  for  the  full  GRPA, 


2-49) 


the  lack  of  fulfillment  of  equations  (2-41,2-42,2-44)  can  be 
attributed  to  the  Y-component  of  Q+  in  equation  (2-49). 
These  components  contribute  to  the  construction  of  the  B- 
block  (off  diagonal  blocks)  of  the  propagator.  The 
generalised  TDA  has  been  applied  to  a range  of  chemical 
systems  (32),  with  promising  results.  The  extent  of 
contribution  due  to  the  Y-component  has  been  tested  on  LiH, 
Lij,  and  Be.  The  results,  which  comprise  the  fourth  chapter 
indicate  that  the  norm  of  the  vector  Y is  virtually  zero  for 
these  systems.  This  is  reassuring  in  establishing  GAGP  as  a 
good  reference  state  for  the 


Polarization  Propagator. 


functional  Fora 

For  simplicity  we  will  concern  ourselves,  in  this 
chapter,  with  reference  states  that  have  a single  unpaired 
electron.  Examples  of  this  type  are  alkali  metal  atoms  such 
as  lithium,  sodium  etc.  and  molecules  such  as  CH,  ArH  etc. 
Certain  doublet  excited  states  of  the  nitrogen  atom  may  also 
qualify  as  a reference  state  for  that  system.  Without  loss 
of  generality  we  will  consider  our  reference  function  to 
have  m - +1/2.  For  clarity  and  convenience  we  will 
represent  the  orbital  which  is  singly  occupied  in  the 
reference  state  by  ♦ . The  open  shell  GAGP  reference  state 
is  of  the  following  form: 

| GAGP>  - (cfiTTkr1  °As<G+)k(*2k+l*2E+I  ••••♦s:i»o>lvac>3 

The  calculation  and  the  optimization  of  the  reference  state 
require  the  same  formalism  as  in  the  closed  shell  case.  The 
only  minor  difference  occurs  in  the  construction  of  the  two 
matrix,  where  the  diagonal  element  in  the  D-block, 
corresponding  to  * is  zero,  and  all  terms  in  the  C-block 
using  ♦-  are  also  zero.  Also  while  using  equation  (2-23)  for 
obtaining  the  gradient  elements  and  equations  (2-29)  to 
obtain  the  hessian  elements,  variations  involving  ♦-  are 


Single  particle  excitations  out  of  the  doublet  reference 
state  may  generate  either  a pure  or  mixed  spin  states  which 
are  the  following: 

1)  A singlet  excitation  operator  generates  a doublet  state, 

2)  a triplet  excitation  operator,  yielding  dm  - 0, 
generates  a mixture  of  doublet  and  quartet  states  of  m - 
1/2, 

3)  a triplet  excitation  operator,  yielding  dm  - 1,  will 
generate  a pure  quartet  state  with  m - 3/2, 

4)  a triplet  excitation  operator,  producing  dm  - -1,  will 
generate  mixed  doublets  and  quartets  of  m - -1/2. 

While  one  could  choose  case  1)  and  3)  to  obtain  pure  spin 
states  an  additional  complication  arises  for  doublet  excited 
states.  A singlet  excitation  operator  will  essentially 
generate  three  unpaired  electrons  out  of  the  open  shell  GAGP 
(3-1),  with  mg  ■ 1/2.  In  such  a case,  the  doublet  manifold 
is  two-dimensional  (45),  and  any  linear  combinations  of  two 
linearly  independent  spin  functions  from  this  space  are 
eigenfunctions  of  S . Either  a branching  diagram  scheme  (45) 
or  Young  Tableaux  diagrams  [46]  may  be  used  to  construct  the 
proper  eigenfunctions.  Possible  functions  that  span  the  spin 


11  *ij,i  * l*o*r*jl  * l+o+i+i I «»-*•} 

21  *i),2  ' *l*o*l*jl  - I *o*I* j I - l*0*3*i I |3-2b> 


For  all  possible  i,j 


(3- 


L yield  <3-2a). 


VtV  ■{>])  - 9i**j#i+  aJaI> 


- 9j aIa5>  - a5aI> 


s - V2),  is  o£  the  £ 


“™*  ■ «.  -i-j  * ».  ■;*! 


(3-4) 


(3-5) 


<GAGP|(ga<2.Jag.#a5  - (,£.,  - .g.;),  - g^agaga,.;;  - 
(aba.  ' a6ai>>]  kh (^c * adac  * aaa5>  - 9d«ai»d  + a|»a)0l<3*GP> 

- 2^  ^ ^((N-l)  1 < hik* j 1 _ hil*jk  ■ hjk*ll  + hjl*ik)  + <i3llkl>)* 

<GAGP|(ga(2a^aaa5  - (a+a,  - aga;))  - 9fc2agagaaa5  - 

«aSaa  - aBa«>>)atajalak(9o'adac  * aa»5>  - Va>d  + alaa> ))  |G*GP> 


- aoaEaaaaa3#aacao I GAGP>  *Bc  } 

+ <5-HBo>(<«-P|.:4.:.eVo-  aoaBacadaaacaBao  lGAGP>*ad^} 

- 9bgc{<ob||ao>(<GAGP|a£a±aja0  + aja^agaga,,  + a*a±a£abaga0 

♦ aoagabaaaaabasao  lGAGP>‘sa 

- <oc||io>(<GAGP|.+.±a+aeaja0-  a*a^a*a^jabac.ga0 1 GAGP> ibd) } 

- 9a9d(f “ I ‘ E0>  (<GAGP I «oaBaaao  ♦ aoaBa£aca3ao  + aoa6'Vaaaao 

- <^MEo>(<GAGP|,;,^adaBa0-  a^a^a.a^ajGAGmJ) 

+ gb9d{<5bllSo>^<GAGplaoa5aaao  * *o*gacacaaao  * aoaiababaaao 

- ao*iabasacabaaao I GAGP>  4sa  ) 

+ <°dMa°>(<GAGP|a^a5ad.5a0-  »Ha^badaiao  l0'''31'^} 
(3-6) 


9a'  «b'  9C. 


diagonalizing  the  hamiltonian  matrix  of  the  above  type  would 
then  in  general  be  a linear  combination  of  the  functions  *, 
and  *2*  *£  the  two  sets  of  functions  did  not  interact 
through  the  hamiltonian,  then  each  set  of  eigenvalues  and 
eigenfunction,  obtained  by  diagonalizing  the  hamiltonian 
would  be  independent  of  the  other  set.  This  would  have  meant 
losing  certain  eigenvalues  by  not  considering  both  sets  but 
not  affect  the  exactness  of  the  results.  Because  of  equation 
(3-4),  eigenvectors  obtained  with  any  one  set  of  operators 
will  not  necessarily  correspond  to  proper  excited  states  of 
the  system. 


GTDA  And  GRPA  For  Doublet  Excited  States 


Obtaining  the  quartet 

excited  stat 

es  is  quite  similar  to 

obtaining  triplet  excited 

states  in 

the  case  of  singlet 

reference  states  and  therefore  will 

not  be  included  in  the 

following  discussions.  To 

carry  out 

a GTDA  calculation  for 

the  doublet  states  we  wil 

1 need  to 

consider  the  following 

the  types  of  integrals. 

1)  <gagp|  [q1E>1>  , [kh  , 0'D 

1 ] ] | GAGP) 

(3-7) 

2)  <GAGp|[q(Di1  ,[kh  , 0<D 

2>  ]]|gagp> 

(3-8) 

3)  <gagp|[q<D2>  ,[kh  , q<D 

1>+]]|gagp> 

(3-9) 

4)  <GAGP|  [Q ( °2  ' , [kh  , Q(° 

(3-10) 

The  excitation  energies  may  be  obtained  by  diagonalizing  the 
entire  matrix,  while  reviewing  the  consistency  of  the  GTDA, 
the  following  problems  are  encountered.  In  order  that  we  may 


calculations,  the  excited  states  generated  by  (3-3a)  and  <3- 
3b)  must  satisfy  the  generalized  Brillouin  condition,  ie 


where  and  Cj  are  expansion  coeffients  corresponding  to 
excited  states.  The  stationarity  condition  imposed  on  the 
GAGP  in  obtaining  the  optimized  reference  state  ensures  that 
it  satisfies  the  generalized  Brillouin  condition  with 
respect  to  the  first  doublet  excitation  operator,  ie 


exclude  the  reference 


<gagp|[kh  , CjO*01  >+  + c2q<I>2>+]]  |gagp>  - 0 


Kg AGP | [kh 


(3-12) 


which  then  leaves  us  to  enquire  if 


<GAGP | [kh 


(3-13) 


Expanding 


(3-13), 


Wlis:  <ij||ki>,;»;«1,k , - .s.s)) 
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“ 4 gb<o6||ao>  <GAGP |a*agaaaaaBaQ-  a*agaBa0 | gagp> 

+ 4 ga<oa||b6>  <GAGp|a*ala*aba5a0-  a Ja Ja5aQ | GAGP> 


(3-14) 


and  |*2)  is  small  and  that  the  eigenvalues  and  eigenvectors 
will  not  be  affected  too  much  by  this  approximation.  In  this 
particular  formalism,  we  need  to  include  in  the  set  of 
excitation  operators  (otjfDj)),  of  equation  (3-3a),  the 
excitation  operator  , which  involves  excitations  of  the 
open  shell  spin  orbital. 

<l-(l-v.)>llj/2)1/2(a}a0  - gja£a3)  lSjSk  (3-15a> 

ajag  2k+l<a<s  (3-15b) 

The  adjoints  of  the  above  operators  annihilate  the  GAGP 
function. 

We  applied  GTDA  to  the  Li  atom,  using  a 5s  2p  GTO  basis 
set  (48]  and  a larger  basis  set,  obtained  by  augmenting  the 
the  basis  set  with  a diffused  d Gaussian  type  function.  The 
results  for  the  first  three  excitation  energies  are 
presented.  We  compare  our  results  with  experimental  findings 
(47),  and  also  other  theoretical  calculations,  using  the 
propagator  approach.  These  include  calculations  by  Swanstr^m 
and  Jorgensen  [19],  who  have  used  50  STO's  to  evaluate  the 
first  seven  excitation  energies.  They  have  used  both  the  HF- 
RPA  ( TDHF ) , and  the  second  order  polarization  propagator 
approach  (SOPPA)  to  obtain  their  results.  In  the  SOPPA 
calculation,  the  two-particle  two-hole  propagator  elements 
required  in  the  generation  of  the  degenerate  set  of  doublet 


results  indicate 


states  (3-2b),  are  incorporated  as 
perturbative  correction  to  the  tdhf.  The 
the  effect  of  this  correction  to  be  quite  small.  The  first 
three  excitation  energies  calculated  by  AGP-TDA  are  quite 
accurate.  However  to  obtain  higher  excitation  energies  one 
must  use  larger  basis  sets.  The  higher  2s  excitation 
energies  can  also  be  further  improved  by  the  inclusion  of 
the  N-particle  diagonal  operators,  equation  (2-41),  which  is 
now  in  the  process  of  being  coded.  Tables  1 and  2 present 
the  basis  sets  used,  and  table  3 contain  the  Li  results. 
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Type 


0.001367 

0.010425 
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0.160701 

0.344604 
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1.00000 


0.03877 

0.236257 


3.010425 

3.049859 
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AGP-TDA,  AGP-RPA  BASED  RESULTS 
AND  CONCLUSIONS 

Potential  energy  surfaces  for  ground  and  excited  states, 
and  other  spectroscopic  properties  such  as  transition 
moments,  and  radiative  lifetimes  for  LiH  and  Lij  have 
already  been  successfully  produced  by  AGP-TDA  (44,49].  There 
is  virtually  complete  agreement  of  the  shapes  of  curves 
produced  by  AGP-TDA  for  LiH  with  extensive  NCSCF 
calculations  (50J.  The  Lij  results  are  also  very  impressive 
and  compare  well  with  experimental  findings  as  well  as  the 
leading  theoretical  calculations  (references  in  (49)]. 
However,  the  strength  of  AGP  as  a reference  function  would 
not  be  fully  understood  without  similar  calculations  using 
AGP-RPA.  The  contributing  terms  of  the  "A”  block  (following 
the  notation  from  chapter  one)  are  the  following 


1)  COlq.jHqj^lOi 

for  all  possible  i,j,k,l  (4-1) 

2)  <o|q1Jq£1*|o> 


The  first  term  is  the  hamiltonian  expectation  value  between 
two  single  particle  excitations  out  of  the  GAGP.  The  second 
term,  by  virtue  of  the  equality  that  was  proved  by  Jensen, 
Weiner  and  Ohrn  (27]  and  Sangfelt,  Kurts,  Elander  and 
Goscinski  [44],  is  a constant  quantity  and  is  equal  to  the 
reference  state  energy  for  the  diagonal  elements  and  sero 
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foe  the  others.  The  AGP-TDA  (which  is  equivalent  to 
evaluating  the  A block  of  the  propagator)  is  therefore 
essentially  a mono-excited  Cl  matrix  which  on 

diagonalization  will  yield  energy  differences  between  ground 
and  excited  states  (which  are  also  eigenstates  of  the 
hamiltonian) . The  contribution  to  the  B block  of  the 
propagator  arises  from  the  following  term  : 


:0lqi  j'JklHl°> 


which  represents  interaction  of  the  reference  state  with  two 
particle  excitations  out  of  the  reference  through  the 
hamiltonian.  The  B block  which  is  included  in  the  RPA  matrix 
constitutes  a measure  of  the  contamination  of  the  reference 
into  the  excited  states.  In  the  following  sections  the 
comparative  results  for  AGP-TDA  and  AGP-RPA  for  LiH,  Li-  and 
Be  atom  are  presented  together  with  similar  calculations 
based  on  a HF  reference  state,  since  the  computer  codes  for 
evaluating  the  diagonal  operators  have  not  been  developed 
yet  for  the  RPA  scheme  we  have  excluded  those  from  TDA 
calculations  as  well  for  proper  comparisons. 


Potential  Energy  Surfaces  (LiH  and  Li-) 

The  low  lying  singlet  and  triplet  excited  states 
together  with  the  reference  are  portrayed  in  figures  1-8. 
For  both  systems,  AGP  based  TDA  and  RPA  do  not  exhibit  any 


detectable  difference,  and  for  each  state  the  curves  are 
completely  superimposed.  The  weakly  bound  (a  well  of  .035 
ev)  characteristic  of  the  B state  of  LiH,  as  predicted  by 
experiments  (51J,  is  not  yielded  by  the  AGP  with  this  basis 
set.  However,  the  HF-RPA,  HF-TDA  predicts  quite  a deep  well 
of  1.5  ev.  While  both  singlet  excited  state  curves  of  LiH  in 
the  HF  picture  seem  to  rise  at  large  internuclear  distances, 
the  shape  of  the  A *E+  is  slightly  improved  by  the  HF-RPA. 
The  2 3E  and  3 3E  of  LiH  are  predicted  to  have  avoided 
crossings  by  both  AGP  and  HP,  however  the  HF-RPA 
calculations  run  into  a serious  stability  problem  for  the 
triplet  states  (discussed  in  chapters  one  and  two)  and  the  1 
3E  state  breaks  down  completely.  An  identical  trend  is 
observed  for  the  Li-  molecule.  The  shape  of  the  A 1E+  state 
is  slightly  improved  as  before  by  the  HF-RPA  over  HF-TDA. 
The  triplet  instabilities  are  just  as  pronounced  for  Li,  as 
in  LiH  and  as  a result  of  this  the  HF-RPA  curve  for  the  b 
\ state  cannot  be  plotted.  Thus,  comparing  AGP  and  HF 
results  for  the  triplet  states  show  that  a remarkable 
improvement  has  been  achieved  for  the  AGP  regarding 
stability.  We  also  present  in  tables  6 and  7,  numerical 
results  indicating  the  effect  of  the  RPA  in  comparison  to 
TDA.  In  these  tables  u denotes  the  excitation  energy  between 
the  ground  and  excited  states. 


0 1 59613 
0.07455 
0.02867 


7042 

56504 

938258 


0.372145 

0.155838 


0i01621 

o!24579 

0.47019 

0.34540 

0.003509 

1.00000 

1.00000 


0.00894 
0 ! 94535 


0.000895 

0.065862 

0.178008 


1.00000 


1175 


57 


Type 


0.001367 

0.010425 

0.049859 

0.160701 

0.344604 


1.00000 

1.00000 

1.00000 


0 1 2667 
0.07201 


1.00000 


Energy  (Har.). 


ss 


Energy  (Har.). 


59 


Energy  (Har.). 


-7.700 


60 


ilatic 


Energy  (Har.). 


61 


Energy  (Har.). 


62 


Energy  (liar.). 


63 


Energy  (Har.). 


Calculi 


Energy  (Har.). 


65 


Figure  8.  Lij  Triplet  States  (RPA,  TDA  Calculation  with  HF) 
+ RFA  Calculation 


As  pointed  out  in  chapter  two,  the  contribution  of  the 

vector  Y,  (2-51).  We  therefore  estimate  the  extent  of 
consistency  of  the  AGP  reference  state  by  1)  the  smallness 
of  the  vector  Y and  the  2)  the  relative  magnitude  of 


matrices  'A'  and  'B'  ie  the  magnitude  of  the  ratio  of  the 
traces  of  the  squares  of  matrices  A and  B.  This  ratio 
indicates  that  the  dominant  (diagonal)  contributions  of  the 
'A'  block  are  several  orders  of  magnitude  higher  than  those 
of  the  ’B' . For  LiH,  this  ratio  is  larger  in  AGP  than  in  HF 
for  all  the  symmetry  blocks.  The  results,  though  tabulated 
only  for  the  equilibrium  geometry  of  the  ground  state,  are 
typical  of  the  entire  surface.  The  norm  of  the  vector  Y is 
virtually  vanishing  for  the  AGP  indicating  in  these  cases 
t reference  function.  Xt 


is  interesting  to  note  that  the  1E*  state  of  Li2  has  the 
ratio  of  the  traces  of  the  squares  of  A and  B,  larger  for  HF 
than  for  AGP,  which  suggests  that  the  relative  contribution 
of  B matrix  is  smaller  in  HF  for  this  symmetry  block. 
However  the  norm  of  Y vector  in  AGP  is  much  smaller,  which 


also  reflected  in  the  agreement  of  the  plotted 


TABLE  6. 


Transition  Moments 


The  lifetime  of  electronic  excitations  is  a particularly 
useful  quantity  since  it  cannot  be  measured  experimentally 
very  easily.  For  the  calculation  of  lifetimes  one  needs  to 
know  the  electronic  transition  moments,  and  these  are  quite 
easily  calculated  in  the  propagator  formalism  (44].  The 
transition  dipole  M0|t  between  the  ground  and  an  excited  state 

M0k  - <0 | r |k>  (4-4) 


where  r is  the  dipole  length  vector.  Equation  (4-4)  may  be 
rewritten  using  the  excitation  operators  as 


5 

i>j 


cjj  <0 1 1 r ofjlIO 


(4-6) 


(4-6)  is  possible  when  the  vacuum  condition  is  satisfied,  we 
may  express  the  dipole  vector  r as 


where  r is  the  dipole  element  between  orbitals  m and  n. 
Substitution  of  (4-7)  into  (4-6)  leads  to  a very  simple 
expression  for  the  transition  moment. 


- u.) (N.  - Nj))*1 


i>l 


We  only  calculate  the  ground-excited  state  dipole  moments 
and  not  excited-excited  state  dipole  moments.  The  excited- 
excited  state  transitions  can  also  be  handled  in  the  AGP-TDA 
formalism,  but  one  might  question  the  corresponding  RPA 
calculations  since  the  excited  states  are  only  well-defined 
when  the  B-matrix  is  small.  The  shapes  of  the  LiH  and  Li, 
transition  moment  curves  are  in  excellent  agreement  with 
other  leading  MCSCF  calculations,  whereas  the  HF  results  are 
quite  inadequate,  even  though  the  HP-RPA  is  somewhat  better 


Transition  Moments  (a.u.). 


RPA  Calculation 


Transition  Moments  (a.u.). 
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• TDA  Calculation 
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• TDA  Calculation 


Calculation 


Transition  Moments  (a.u.). 


Figure  12.  Li.,  Transition  Moments  (RPA,  TDA  Calculation 
• TDA  Calculation 


Calculation 


Results 


Using  a 6s  3p  2d  basis  (table  8),  we  compare  the  first 
two  excitation  energies  for  S and  P symmetries  for  both 
singlets  and  triplets.  The  agreement  of  the  AGP  results  with 
experimental  values  [ 52 ] is  rather  good  and  there  are  no 
instabilities.  The  results  for  *S  are  in  fact  much  improved 
with  the  inclusion  of  diagonal  operators,  (2-41)  and  yield 
an  excitation  energies  of  6.91  and  8.099  ev.  The  agreement 
of  the  AGP-TDA  and  AGP-RPA  is  also  reflected  by  the 
dimunitive  norm  of  Y and  large  ratio  of  the  traces  of 
squares  of  A and  B.  The  AGP  calculation  results  for 
transition  probability  compares  very  well  with  experiment 
153). 
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APPENDIX  1 • 

GRADIENT,  HESSIAN  AND  PROPAGATOR  ELEMENTS 

The  different  types  of  gradient  elements  areV 
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for  2s+l<k,l<2s+(N-2k) . 


Habcd  ■ u(.»i_«a)Moi.;.;.J.b  • 

, <m.d|bdi  , (i*|b«m<oK.;Vk  * «-a*.  - 

* ab‘Baaad  ' asaaaoac  - aaaiaaad  - aba§aBa«l0>  - 

- £ (i(,b|ndfe,;,i.E.b  . aXaaaa  - aXaban  - aXada. 

• "Kv,  - IaMaaa,l">l*.b 

a I<  an  |nd ) I <0  |aaaEajaa  a anaHaaad  “ anaaaaan  " anadadan 

a XaXab  ' XaEaHan  1 0>l  ‘bn 

a 2(nb|cnK<0|.aaiaBab  • .K-ji,  - aXaban  - aXa«a« 

* 2aaananaa  - 2aaaiananl0>,Sad 
+ 2 ( an  | cn ) ( <0  |a*.ia-«a  ♦ a+afoa,.  - a^a^n  - <a£aean 
+ 2abananab  " 2abaBa5anl0> » sbd 

* 4(ac|nn)(<0|a*a*aaan  + a*a*acan  - 2a+a£abaj0>)«bd 
+ 4(ad |nn) «0 |a*a*aaan  ♦ a*a*ada„  - 2aJ;aJaban|<»)Sbc 
+ 4(bc|nn)(<0|a+ajaban  + «„a*acan  - 2a*a*aaan|0>)6ad 

* .(bd|.nmO|,Xabab  * aXada,  - “K*.a. I "»*..} 


fi  - <Vd  ♦ ; t2 


11  Pabcd  ■ 4(.b|cd.((f1  * fj)<0|a^.d.b+  - «&.c«b 

- aaadadaa I ®>  * 

+ 2(aclbd><£1<0|»bao«cab  - a£a2adab  ' aaacacaa  * aaadadaa 
+ aaa5acac  + abaBaaadl0>  - £2<°Kaiaaad  + aSaBaaad0>> 

+ 2(ad|bc)(£2<0|a*a*acab  - a£a*adab  - a^acaa  + a^a*adaa 

* aaaaa3ad  + aSaBa5ac I °>  ' £i<° la:aiacac  + a£aBaaadl°» 

- £ 

- (u.i»bi«!<»i.;.|.8.n  - 

- £1<® labaBanan * + 2(ad|nn)f2<0|a*«*a„a,-  aja^a,,  |0M  6bc 

- (lnb|cn)(t2<0|a^aK,n  - a^Vb  * aXa„aa|oV 


- *i<® l®g®fanan l0> 1 + 2(bc|nn)f2<0|a*a*anab-  aaa*anaa|0>)Sad 

* ((®nlen)(f1<0|a*.iagan  - a^a,  ♦ »Xabanl0> 

- £2<0.lab®5an*nl0>>  * 2<®c|nn)fl<0 ®n»baban l0>) *bd} 
2>  pibcd  - 

2(ac|bd)(£1<0|a*a+acab  - a£a*adab  - a+a*acaa  + »a»d*d®a 
+ »M«c«c  * abaBaaadl°>  - f2<0|a:®i®aad  ♦ ®baBacacl0>l 
+ >(«ibe)i'j<oK®;*e®b  - ®s®3ad®b  - ®:a:®c®a  * 

* aMaaa<l  * aSaBacacl°>  - £l<°l®iaiacac  ♦ abaBaaadl0>> 

- Y.  j<  (nblnd> I £1<° l*baBan®n  - abananab  + aXanaal0> 

“ £2<0laaa5®5anl°>>  + a(bdh®>«x<0  aXv>a  l°>  > ‘ac 

- ((an|nd)(f2<0|aJai.Han  - aja^a,  + «Xab*nl0> 

- £l<® l*b®Banan|0>f  + a(-|n.|f2<0|.;<Va-  ®Xabanl°»Sb= 

- ((nb|cn)(£2<0|aJagaHan  - a£a*®nab  + a*a*»„aa|0> 

- £l<0|aM®Han|0»  ♦ 2(bc|nn)£2<0|a^«nab-  aXan®a|0»»ad 

* * ta»i I cn)  (£i<0|aaaj«g«n  - »Xa„aa  + aXabanl0> 

- £2<0labaBananl0>)  + 2<ac|nn)£1<0|a*a;anaa-  a+ajaban | 0> ) «bd) 


«1)  <0|tgb.;»g«aa5  * 9t»>Jab»;  [K„  , 9da|a+a5a0 

+ 9ca5adacao  H|o> 

.2,  <0|Cgb(.Jaa  - a|a-J  + gja^  - at«B]IV  ,„(«>,  - .±.,) 
* 9c(adac  *aac)ll|o> 

a3)  <0I  '9baoaBaaao  + ®aaoaaabao  >kh  • 9d(a>d  ~ a5aa) 

+ 9c  (adac  “ aaac)"|o> 

°4 ' <0|'9b(abaa  - aBaa)  + Safa^h a5aB)>KH  < 9da5adaEao 
+ 9caoada5ao  'll<» 

.5)  <0|fgbaj;agaaa5  - g^a^g  (kh  , gd(a*«d  * a|ag) 

- 9c(adac  + aaas))i|o> 


<0H9baJa^aBao  ♦ gaaia*as«Q  (kh  . 9da±a+aja0 
+ 9caoadacao  > » |o> 

b2)  <0|l9b(aaab  - aSaB)  + 9a (abaa  " 9d(a^  ' a5aa) 

+ 9c(a3ao  - aa*c)|1l°> 


85 

-3.  <0 1 1 9baoaaaEa0  ♦ 9aaia*a5*o  1«h  ' *d(«^  “ aSaa) 

+ ®e (®d*o  - »aae)>ll°> 

b4>  <0|,9b(aaab  * *iaB)  ♦ 9,Kaa  - aBas)>KH  ' ®daoadacao 
+ ®ca5adacao  m°» 

b5)  <0|lgbaia*aBao  + 9aaia*a;ao  [kh  , gd(a*a„  + •§.,) 

- 9c(adac  ♦ •3ae),,l«> 

The  expression  for  each  term  is  as  follows 


(ac |bd) (gbgd  <0 |a*a£a+a±a5a8aBa0 

' aoaEabaaaaada3ao 
- aoaBacacaBaol0> 
®a9cl<°laoaiadaaaBabaIao 


- 9a9d(<°laoa5a^aBaBabaaao  " aoaiaaababadaaaol0>> 

" ®b9c«°laoaBadaaaaaaaBao  " aoaBabaaaaacacao l0> ») 
+ <adlbc>(<9b9d  - 9a9c)«0|a*a|a*a|a;iaaBao 


»aaol°> 


* 3agd«0|a*a|a;acas.o|0>.  ♦ 9b9c«0  la^a^  | o> 

+ (9a9d-9bgc)«0|a;aia;a|a5abasao  ♦ vKa£abacaSao 
+ aoai*baaaaaba5ao  - aoa5aaababadaaaol0>>) 

+ (,°leo,(9b9d<0laoaKacaBaol0> 

*<3b9d  - 9a9c * f <® I aoaBaaaeacaaaBao  + aoaBa£a>aaaaBao l«» > *fia 

+ ,boldo,Kv0laoaSacaeaaaol°> 

+(9*9=  - 9b9d)«0|a*aga*a|agaaaEao  * •^gaJ.X^.g.Jem,. 

+ ( 80 1 do ) (9a9d<° I aoaBadadaBao I 9> 

+ <9a9d  - 9b9e)(<°laoaEa^adadaaaBao;+  aoaBaaadaaaaaBaol0>,Sbc 
+ (bo|co)(gbgc<0|a*aia*aeasao|0> 

*(9b9e  - 9a9d)«0|a*a±aSa*aeaba5ao  * a^a^ajaJO)^ 

+ Z<an|cn,(-  9b9d<°laoaBaXaBao|0>  + 9a9c<0Kaaa^anaaaol0> 
+(9b9d  - 9a9c * * <°  I aoaBaaananaaaBao  " aoaBaKaaa5aBaol0>»‘bd) 

bn|dn)(-  9a9c<0|a*a^ana5a<)|0>  + 9b9d<0 la^^a^g^ [o> 
- 9b9d><<°Ka6aXanaaaBao  - aoaBabaaaaanaHaol°>»‘ac) 

+ J3(»n|dn)(-  9b9c<0 1 aoaBananaBao 1 0>  + 9a9d<° l<aKanaaao l0> 


♦(9b9c  - 9a9d,(<°laoaiaaaIanaaaBao  " aXaXaaaiaBaol0>1‘bc) 

+ Dbn|en»(-  9«9d<°laMananaSaol0>  * ?b9e<0 |a^BananaBao l0> 
+(9a9d  - 9b9cl I <0 l*o®BaaananaaaEa0  * aXabaaaaana5aol0>  1 ‘ad) 

+ (ac|oo)(gbgd  <0|»^aB«o|0>  - gage  <0|a*a±ajao|0> 

-<gbgd  - 9,9C ) <0 | »o*Baaaa*5ao I °>) ‘bd 

+ <bdl°°)(-9b9a  <°laoaBaBaol0>  + 9a9c  <°laoa;aiaol0> 

♦ '9b9d  " 9a9c)<#la0aBaaaaaBaol0>)‘ac 

+ (*<3|oo)  (gb9c  <0|a^a5aJO>  - gagd  <0  la^a,,  1 0> 

-(9b9c  - 9a9d)<0|.;.^.aaB.o|0>)Sbo 

* (bc|ool(-gb,d  <0 |a*aBaB»0 | 0>  + g#gc  <0|.;«±,5,o|0> 

-(gbgd  - 9,9cX0|.;a^aaaBao|0>)6ad 

*clnn,(  9bV0!*^KanaBao!0>  - 9a9c<0 1 aoaSanana5ao I 
- 9a9c ) <«> I aoaBaaananaaaEao l0>) ‘bd) 

♦2^<bd|r,n>(-  9b9d<0|a+»ga;.n»Baol0>  + 9a9e<°  laoaKanaaao  l0> 
+ (9b9d  - 9a9c>«°laoaBaXanaaaBaol°>»‘»c) 

+2^(ad|nn)(  9b9c«>  la^a^a,^  |0>  + 9a9d«>  l«XaXaaaol0> 


-<9b9c  - 9a3d ) < <0 1 aoaEaaananaaaEao I °> > *bc) 

+aZ(bc|nn,("  ^b®c<0 laoaSan*naBao l0>  + «.V° laoaKanaiaol°> 
*(9b9c  - 9a9d) ( laoaBaaananaaaBao |0> * sad) 


(ao|co)(9b9d  <0|a^aBao|0>  +.gagc  <0 |a>|a5a0 |0> 

- 9b9d  <“l*oMaaaaBaol0>  + Wc<0  laoa5ababa5ao  l0>)  sbd 
+ (bo|do)(gbgd  <0|a*agaBao|0>  ♦ ga9c  <0 |a*a|as.0 |0> 

- 9b9d  <0 1 aoaBaaaaaBao  1 0>  + *a  V°  H 4ac 

+ (ao|do)(gbgc  <0  la^aJO^g  * gagd  <0  |.>|«aa„  |0>*bc 

- 9b9c  < 0 I aoaBaaaaaB*o  I ® > *Bg  + W°  la>Kabaaaol°>4bc) 
+ <b°l=°>(9b9c  <0KaBa6aol°>Sac  * 9a9d  «>  |a*a±a5ao|0>*sa 

- Vc  <°laMaaaaaBaol°>Sac  + *aV0  K^bVo  l°>4Sa) 


a4)  - 


( ao | co ) (gb9d « 0 1 adababa0  - aoaaaa*o|0> 

+<9b9d  + 9a9c*<® laoaaaeac*a*o  “ aXacacabaol0>)*bd 
* (*>o|do)(9age«0|a*aJ;aaao  - a*a*abaol0> 

*-9b9d  + 9a9c,<0laXadadabao  " aoaaadadaaaol0>)‘ac 

+ (ao|do)(gb9c(<0|aXabao  ' aSa2a,a0l*> 

+<9b9c  + 9a9d)<0laoaaadadaaao  ' <abadadabao l0>) ‘be 
+ <bo|co)(gagd«0|a*a*aaao  - a+a£abao|0> 

+<9b9c  * 9a9d*<0 laoabacacabao  ' aoaaaJacaaaOl0>Kd 

Tern  a5)  - 

(ao|coj(-gage<0|a;.±as.o|0>  - W° laHaBao l°> 

+ *9a9c  * 9b9dK0|»XaSaba.aol0>)4bd 
+ <a°l«Jo>  ( 9,9d<0laSaiaiaol0>  + 9b9c<0 1 ao*B*B*o 1 9> 

* (9b9c  - 9,9d*<0 laoaaababaaao l9>) *bc 
+ ( bo  I co ) (-gagd<0 1 a^afaja,, 1 0> . - gbgc<0 1 a^aBa„ 1 0> 

+ <9b9c  ' 9a9d»<0laoaaababaaaol0>Kd 
+ <bo|do)(  9a9c<0 1 ao«aaaa0 1 ®>  * 9b9d<0|a^.5.o|0> 


+ «Vc  ' 9b9d><°laoat«b,,b“»“ol0>)‘ac 


(bo|do)(9a9d<0|a*ala5aa  - ♦ adaaaaad  " ^adadaJ0> 

+ 9b®c<® labaBaBab-  aaabab*a  +:a£aSacac  ~ aJadadac l0> > *ac 
+ ( bo  1 00  H 9a9c<0 1 aaaSaaaa  ‘ aaaSabaa  + aMacac  “ aJaSadaol0> 
+ 9b9d<0|a*a|aBab-  a+ajabaa  ♦ ad»aaaad  " <a5adac l0> >*ad 
+ (*° l^°) (9b9d<0 |abagagab  - aja^a,  + adaaaaad  ' acadadacl0> 

+ Vc<0laaaSaaV  aXabaa  + aMacac  “ a^adac  l0>  > Sbc 

+ ( ao | co) ( 9a9d<° laaa5ajaa  - aja^a,  ♦ a^agad  - a+a+adac|0> 
+ 9b9c<0|a*agaBab-  »Xabaa  * a*a±aBac  - a^+adac  1 0>  )S„d 


ao  - a:a:aaaoi°> 


(ao|co)(9a9d(<0|a;a;aba 


*‘*a«d  + 9b9e) <0 | ®o*a*cac*aao  “ aoa£acacabaol0>)Sbd 

* (b°|do)(9b9c(<0|a^aaao  - «^S|0» 

+(9a9d  ♦ 9b9c)<0|aJa‘a*adabao  - a*aX*daaa0 1 0>)  iae 

♦ (ao|do)(gage(<0|a*a*abao  - aJa^aJO) 

+<9a9c  * 9b9dX0|a^,*adaaao  - a^a+ad»bao|0>)*bc 

+ <bo|co)(gbgd«0|»^.aao  - aXabao|0> 

+,»a«c  + ®b9d X 0 1 ao*bacacabao  " aoaaaXaaaol0>)*ad 

Tern.  bS,  - * 

(bo | do) (gagd<0 |a*a|aj.a  - .;.Jab.,|0» 

- 9b9c<0 I abaBaB*b”  aaababaal°>,5ac 
+ <bo|”>'-9aV°Kaiaaaa  ' aXabaal°> 

+ 9bV°labaBaBab-  aaababaal°>,4ad 
+ (ao|do)(gbgd<0|aJ«|»Bab  - a+a£abaa  1 0>- 

- V>c<0laaaiaSaa-  aXaba»l°>)5bc 
+ (ao|eo)(gagd<0|aaaja5aa  - aa«£abaal0> 

- W*l*Hw  aaababaal0a  )Sbd 


BIBLIOGRAPHY 


1.  J.  Linderberg  & N.  Y.  Ohrn,  "Propagators  in  Quantum 
Chemistry"  (Academic  Press,  London,  1973) 

2.  R.  P.  Feynman,  Rev.  Mod.  Phys.  ^0,  367(1948); 

J.  Schwinger,  Proc.  Nat.  Acad.  Sci.  U.S.  37,  452(1951) 

3.  T.  Matsubara,  Prog.  Theor.  Phys.  23,  629(1955) 

4.  D.  N.  Zubarev,  Soviet  Physics  Uspekhi,  Vol.  3(9); 
320(1960) 

5.  P.  M.  Dirac,  Proc.  Cambridge  Phil.  Soc.  26,  376[5.10| 
(1930) 

6.  D.  Bohm  and  D.  Pines,  Phys.  Rev.  92,  609(1953) 

7.  P.  Ring  and  P.  Schuck,  "The  Nuclear  Many-Body  Problem" 
(Springer-Verlag,  New  York,  1980) 

8.  D.  J.  Thouless,  "The  Quantum  Mechanics  Of  Many-Body 
Systems"  (Academic  Press,  New  York,  1961) 

9.  E.  Dalgaard;  Int.  Journal  Quantum  Chemistry,  15, 

169(1979)  — 


10.  D.  J.  Rowe,  Rev.  Mod.  Phys.  40,  153(1968);  T.  Shibuya 
and  V.  Mckoy,  Phys.  Rev.  A 2,  2208(1970) 


11.  J.  Linderberg  and  M.  Ratner,  Chem.  Phys.  Lett.  6,  37 
(1970);  J.  Linderberg,  P.  Jorgensen,  J.  Oddershede, 
and  M.  Ratner,  J.  Chem.  Phys.  56,  6213(1972); 

P.  Jorgensen  and  J.  Oddershede,  J.  Chem. 

Phys. 57  277  (1972);  P.  Jorgensen,  J.  Oddershede  and 
M.  Ratner,  J.  Chem.  Phys.  61  710(1974). 


500(1973);' 
307(1977) . 


Shibuya  and  V.  Mckoy,  J. 


13.  J.  Oddershede,  Adv.  in  Quantum  Chem.  11,  275(1978) 

14.  J.  Linderberg  and  N.  Y.  Ohrn,  Int.  J.  Quant.  Chem.  12, 
161(1977) 

15.  N.  Y.  Ohrn  and  J.  Linderberg,  Int.  J.  Quant.  Chem.  15 

343(1979)  — 


16.  J.  Bardeen,  L.  N.  Cooper  and  J.  R.  Schrieffer,  Phys. 

Rev.  108,  1175(1957). 

. J.  M.  Blatt,  "Theory  of  Superconducting  Systems"  in  Pure 
and  applied  physics  12,  (Academic  Press,  New  York,  1964) 


iusrss&sv-  “■*•*”*• -■ 

StmSjST"  -1  ■■  »”■  1- 

SisuBl!  °'  °°“*“***  ,“t'  J'  ““*•  — • 

P.  O.  Loudin,  Phys.  Kav.  a,  139,  357(1963) 

P.  O.  Lovdin,  Int.  J.  Quant,  chan,  4,  331(1971) 

'•  “■  '■  a. 

v.:  Si,rsshwr“-  ■”*  ’•  j- 


ii,;s?sr"  •“ J-  >«■  «-•  < 

5;.K:  '■ 

j;.;;i!i.r)“a  >“•  a. 

{iuiinir"'  J'  ‘“id.  6, 

SiafTHS;;'-  “•  am. 


“■  *“  ■•  ’• ol1”'  ”i- j-  «“*• 

"■  »•  «• 

"•  •■ 

“■  fcjrs»i'w*  ■*■ 

"•  Ki.;^ss^.^jr«ri.r&7 

“^^SSriS®S!KS-* 

“•  '■  J-  s. 

”•  Sal;  K'fe^nJnr.s  *— 

**•  *•  “”pt  •""  -■ »■  *•  ■*»•  * 2. 


BIOGRAPHICAL  SKETCH 


Rina  Roychowdhury  was  born  and  raised  in  India.  She 
graduated  from  Ashok  Hall  Girls'  High  School,  Calcutta, 
obtaining  her  high  school  diploma  in  science  from  the  West 
Bengal  Board  Of  Secondary  Education.  She  obtained  the 
Bachelor  of  Science  degree  with  honors  in  chemistry  from  the 
University  of  Calcutta,  she  is  married  to  Shantanu  Basu  and 
has  one  son,  Shourya. 


95 


rmmsmzmm* 


si 


